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^p" 1 ' Abstract. We study the class of open continuous-time mechanical particle systems 

, introduced in the paper by Khanin and Yarmola [llj . Using the discrete-time results 

C$ ' from |llj we demonstrate rigorously that, in continuous time, a unique steady state 

exists and is sub-exponentially mixing. Moreover, all initial distributions converge to 
the steady state and, for a large class of initial distributions, convergence to the steady 
state is sub-exponential. The main obstacle to exponential convergence is the existence 
of slow particles in the system. 

^t" . 

AMS classification scheme numbers: 82C05, 60J25, 37H99 

O 

■ 1. Introduction 

Non-equilibrium steady states for open mechanical particle systems have been 
^ extensively studied in recent years as a part of the efforts to derive the Fourier Law 

from the laws of microscopic dynamics [U El [El [HI [121 [El EE] • The existence of 
non-equilibrium steady states has been shown for only a handful of systems using either 
spectral-theoretic arguments [H El [8] or Harris' ergodic theorem [HJH5]; these methods 
automatically yield exponential mixing and exponential rates of convergence of initial 
distributions to the unique steady state. For systems that do not mix exponentially fast, 
which is frequently the case for continuous-time mechanical particle systems driven by 
heat baths, the methods to prove the existence of non-equilibrium steady states are 
virtually nonexistent. This paper, together with [TT], provides a rigorous proof of the 
existence of a unique steady state with sub-exponential mixing rates for a class of random 
billiards driven by thermostats. 

In the system we consider, particles move freely between collisions with disk-shaped 
random scatterers arranged periodically on a torus. Upon the collision of a particle 
with a scatterer, an energy exchange occurs, in which the tangential component of the 
particle's velocity is randomized according to a Gaussian distribution with variance 
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proportional to the temperature assigned to the random scatterer, and the normal 
component of the velocity changes sign. 

The disk-shaped random scatterers act as thermostats in the following sense: 
Assume we set all the random scatterers at temperature T and start with an ensemble 
of N particles at temperature T 7^ T, i.e., the probability that a given particle has 
kinetic energy within dE near E is approximately ce~^° E dE, which corresponds to the 
Gibbs distribution at temperature T = We show, inter alia, that the distribution 
of the particle kinetic energies converges to the Gibbs distribution at temperature T. 

The particles in this system do not interact, so the analysis can be simplified by 
considering a system with only one particle. In this situation, it is convenient to study 
the discrete dynamics of particle collisions with the boundaries of the thermostats. Upon 
collisions with thermostats, the normal component of the particle's velocity represents 
the memory kept in the system, while the tangential component gets fully randomized. 
Therefore, one can define a discrete Markov chain in position and normal velocity 
variables that represents the original dynamics at collisions. For this Markov chain, 
the existence and uniqueness of the invariant probability measure, which describes the 
steady state, along with exponential rates of convergence of the initial distributions to 
the invariant measure were obtained in [11] via Harris' ergodic theorem. 

In this paper we extend the discrete-time results of [llj to continuous time. In 
section [2] we introduce the model and state the results. We present the proof of the 
main theorem in section [31 Harris' ergodic theorem does not directly apply to the 
continuous time system, since the system does not mix exponentially fast. However, 
one of the key observations in our proof is that the potential (Lyapunov) function and 
the minorization condition required for application of Harris' ergodic theorem in the 
discrete case can be used to obtain upper and lower bounds on certain integrals and the 
measures of certain sets. 

Before proceeding with the actual proof we briefly review the results in [llj in l3.ll In 
I3.2l we prove the existence of a unique invariant probability measure (the steady state) by 
constructing a suspension flow over the discrete dynamics and showing that the average 
of the roof function is finite. The mixing and convergence of initial distributions to the 
invariant measure are obtained by applying [T4l Theorem 6.1] by Meyn and Tweedie for 
general state Markov chains. The conditions of the theorem are ensured by adapting 
some of the arguments in [11] to continuous time setting. We outline the modifications 
required and we refer the reader to [H] for more details. In this sense this paper serves 
as a companion paper to |llj . 

The main reason for sub-exponential rates of mixing is that, from time to time, the 
particle acquires low speed, which leads to long traveling times between collisions with 
thermostats. In 13.41 we show that the (invariant) measure of the set of particles that 
would not collide with a thermostat in time r is of the order of 1/t 2 . This implies that 
mixing rates are, at best, polynomial. 
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Figure 1. Geometric Configuration 
2. Settings and Results 

Let r = T 2 \ U p i=1 Di be a bounded horizon billiard table with a finite number of non- 
intersecting circular obstacles D ly ■ ■ ■ , D p with radii Ri, • • • , R p . Figfj] illustrates such 
a configuration. Non-interacting particles move freely in T with constant velocities 
between collisions with the obstacles. Each obstacle Di, 1 < i < p, is assigned a certain 
parameter which represents the inverse of its temperature. The obstacles play the 
role of thermostats. 

Upon the collision of a particle with a thermostat D iy let v~ = (t>^,t> t _ ) be the 
decomposition of the particle's velocity into its normal and tangential components 
with respect to the boundary dDi. After the collision, the normal component of the 
velocity changes sign, so that v\ = — fj, while the tangential component is absorbed 
by the thermostat and a new tangential component vf is drawn from the distribution 

e~^ lV ^dvf. The outgoing velocity has the decomposition v + = (v±,vf). 
Since particles do not interact with each other, we can reduce our discussion to 
studying a system of only one particle. The phase space of such a system is 

fl = {z = (i,v):i6r,!)GR 2 \{0}}/~ 

where ~ is the identification of points on the collision manifold that correspond to 
incidence and reflection with the choice of the incoming velocities v~ = (v^_,v^). Note 
that if the particle originates with v ^ 0, the probability to acquire zero velocity is 0. 

The dynamics of the system is thus described by a Markov Process $j on Q, which 
is deterministic between collisions with thermostats with random kicks at the moments 
of collisions with dT. Let V 1 denote the transition probability kernel for $ t . 

Definition 2.1 A measure on VL is called invariant under $ 4 if for any Borel set 
Actt and t>0 

fi(A) = (VMA) := / V t (z,A)d^(z). 
Jn 
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Definition 2.2 The Markov process $ t is called mixing if an invariant probability fi 
exists and if, for all Borel A and B, 



J B 

Definition 2.3 An initial distribution v is said to converge to the invariant probability 
measure fi if 

lim ||PV - = lim sup \V l v{A) - fi(A)\ = 0, 

and does so exponentially fast if there exist C, a > such that \V l v — < Ce~ at . The 
norm || ■ || is called the bounded variation norm. 

The main result of this paper is as follows. 

Theorem 1 There exists a unique invariant probability measure for the continuous 
Markov process $ f , which is absolutely continuous with respect to Lebesgue measure. 
This measure is mixing, but not exponentially mixing. All initial probability distributions 
converge to the invariant measure, while convergence rates are at most polynomial for a 
large class of initial distributions. 

3. Proof of Theorem CD 

As a preparation for the proof, we outline the settings and the results of [TT] in I3.ll In 
addition we mention several important ideas from the proof of the main theorem in [TT] 
and state two propositions, to be used later, whose proofs are natural modifications of 
the proofs in [TT]. In l3.2l we prove existence of an absolutely continuous (w.r.t. Lebesgue) 
invariant probability measure for the Markov process $ f . The uniqueness is guaranteed 
by absolute continuity and one of the aforementioned propositions. In 13.31 we show that 
the invariant measure is mixing and all the initial distributions converge to it. In 13.41 
we demonstrate a class of initial distributions for which the rates of convergence to the 
invariant measure are sub-exponential. 

3.1. Review of settings and results of [77] / 

To study the continuous-time Markov process $f it is convenient to consider a discrete in 
time dynamical system of particle collisions with the boundary dT of V. When specifying 
$ f at collision points, we chose to keep old, incoming velocities. This choice enables a 
2-dimensional description of the discrete dynamics characterized by the collision points 
r n e dT (parametrized by arc-length) and the absolute values of the normal velocity 
v±(n) at the moment of collision. Indeed, using a randomly generated tangential velocity 




We say that [i is mixing exponentially mixing if there exist C, a > such that 
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Vt(n) drawn from the probability distribution y ^e~ /3iVt ^ n ' >2 dvt(n), one determines the 
outward velocity and hence, the next point of collision r n+1 and the next normal velocity 
v±(n + 1). This procedure defines a Markov chain <3> on Q = dT x [0, oo); denote it's 
transition probability kernel by V. 
The main result of [H] is: 

Theorem 2 The Markov Chain $ admits a unique invariant measure \i, which is 
absolutely continuous with respect to Lebesgue measure. Furthermore, there exists a 
non-negative function V on Q, as well as constants C > and 7 G (0, 1) such that 

sup ^((r,^),^) -fi(A)\ < C 7 n (l + V(r,v ± )) 
Acn 

This theorem has an important corollary: 

Corollary 3.1 The discrete-time Markov chain $ is exponentially mixing and if a 
probability measure A on Q satisfies + V(r, v±))d\ < oo, then V™\ converges to \x 
exponentially fast with the following bound on the rate of convergence: 

\\V:\-fi\\<CT [ {l+V(r,v ± ))dX. 

The function V above is called a potential or a Lyapunov function and provides the 
rates with which the dynamics converges to the 'center', represented by a level set of V. 
In [TT] V is given explicitly and happens to depend only on v±. A useful observation for 
our analysis is that, by [9, Theorem 3.2], f V(v±)dfi < 00, from which we obtain upper 
bounds on the density of \x near and 00. 

Another important property of the discrete-time system, and the second condition 
for the application of Harris' ergodic theorem, is the minorization condition on compact 
sets. It guarantees the coupling required for fast mixing. For any < v ™ in < t>™ ax < 00, 
let C := {(r, v±) : v ™ in < v j_ < v ™ ax } and let mc be the uniform probability measure on 
C. It was shown in [TT] that there exist iV and rj such that 

inf V N ((r,v ± ),-)>7imc. (1) 

(r,v ± )eC 

The minorization condition for discrete-time dynamics is used in this paper to provide 
lower bounds on the measures of certain sets. 

The proof of the minorization condition in [TT] is accomplished in three steps: (i) 
Prop. 3, which guarantees existence of a regular iV-step sample path from any (r, v±) G C 
to any (r',v' ± ) G C, where C = {(r, v±) '■ v± m < v± < {;™ ax } for some v± m < u™ m and 
-max > yinax. ^ p r0 p ^ which gives bounds on the densities when pushed forward 
along regular paths; and (iii) Prop. 5, which guarantees that for any (r, v±) G C, the 
absolutely continuous component of the two-step push forward of 5( r , Vj _) is nonempty 
and its support contains a neighborhood of (r",i/[) G C. So, originating at (r,v±) G C, 
S(r, v± ) acquires density around (r",v'j_) G C by 5 and the density can be pushed forward 
along the sample path from Prop. 3; and Prop. 4 ensures that we get a uniform lower 
bound on the density at the final point (r', v' ± ). 
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A discrete-time Markov chain $ is called irreducible if for all z G Q, whenever 
Leb(v4) > 0, there exists some n > 0, possibly dependent on both z and A, such that 
V n (z, A) > 0; the definition for continuous-time Markov processes is similar, with n 
replaced by t. One may prove of irreducibility in a manner similar to the proof of the 
minorization condition in [TT] except one needs to deal with sample paths on the entire 
phase space instead of just a compact set C, but there is no need to keep track of the 
lower bounds. The proofs in [TTj of Props 3, 4 and 3 can be modified to obtain the 
following statements about the continuous-time Markov process 

Proposition 3.1 For any two initial states (x,v) and (x',v') in Q, there exists a 
continuous time sample path from (x,v) to (x',v') 

Proof of Prop. 15711 

Let (r, v±) and (r',v' ± ) be the thermostat positions and normal velocities at the 
collision after (x, v) and at the collision before (x', v'), respectively. Then Prop. 3 in [TT] 
guarantees that there exists a sample path between (r,v±) and (r',v' ± ) provided that 
v ±i v '± ^ [ v ±"\ "i"]' However, the only use of the restriction v±,v' ± G [t>™ m , u™ ax ] in the 
proof of Prop. 3 in [TT] is for the bound on the number of collisions along such a path. 
Since the number of collisions is irrelevant for showing the existence of a connecting 
sample path, we conclude that such a path exists for all v± and v' ± (see Lemmas 5 and 
6 in [TT] for details) and this trivially extends to a path between (x,v) and (x',v') in 
continuous time. □ 

Proposition 3.2 The Markov chain $ t is irreducible. 
Proof of Prop. [3721 

Given z = (x,v) G Q and a set A with Leb(A) > 0, let (x',v') be a Lebesgue 
density point of A. Let (r,v±) and (r',v' ± ) be the thermostat positions and normal 
velocities at the collision after (x,v) and at the collision before (x',v'), respectively, 
and let U be an open ball around (r',v'j_). The Props. 4 and 5 in [TT] imply that if 
V -L, v '± £ [f ™ m , f™ ax ], then there exists n such that V n ((r, v±), U) > 0. We are only 
interested in the positivity of V n ((r, v±), U) and do not require any lower bounds, so 
we can drop the assumptions on v±,v' ± . Since the paths from (x,v) to (r,v±) and 
from (r',f^_) to (x',v') are deterministic, we conclude that there exists t > such that 
V\(x,v),A) > 0. □ 

3.2. Existence 

The result of [UJ guarantees existence of a unique invariant probability measure for 
the discrete Markov chain $, which is absolutely continuous. Since dynamics of $ t is 
deterministic between collisions, an invariant measure for $ t may be constructed using 
a suspension flow over the discrete-time chain. This alone does not guarantee that this 
invariant measure would be unique or even finite. We show finiteness by estimating the 
integral of the roof function over the phase space. We obtain an upper bound using the 
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potential V in [11]. Uniqueness follows from the irreducibility of $f, Prop. 13.21 The 
argument for the uniqueness is the same as in [TJ [16] and can be found in the outlines 
of the proofs of the main theorems in either of the two papers. 

For some of our future analysis it is more convenient to use a geometric description 
similar to that used in billiards. Namely, instead of the random variable Vt(n), we shall 
consider a random variable (p(n), which represents an angle relative to the normal vector 
at the point from which a particle leaves after collision. It is easy to show that <p(n) 

is drawn from the distribution p VJ _^((p n )d(p n = \J~^ c l^^ ^ e~ ft ^^ tan2 ^dip n , which 
depends on the fa associated with dDi 3 r and v±(n). We shall also use notation tp' for 
the incoming angle at the point of collision r n+ i. Obviously, r n+ i, v±(n + 1), and ip' n are 
completely determined by r n , v±(n), and (p n . Note that v±(n + 1) = ^|^"| fj_(n). 

Let f(r,v±) be the density of the unique absolutely continuous invariant measure 
for the discrete Markov chain <3>, i.e. dp = f(r,v±)drdv±. Before constructing the 
suspension flow, we need to re-introduce a third variable ip. The invariant measure for 
the discrete process with (r, v±, <p) variables is 

dpx = f(r,v ± ) Pv± ( V )drdv ± d V = f( r ,v ± )J^-^e-^ tan2 ^drdv ± d^. 

V vr cos A (<p) 

Let v by the invariant measure for $ t obtained by constructing a suspension 
flow over the discrete dynamics in a standard way. See [21 p. 31] for details and 
additional references. The suspension flow is defined on a 4-dimensional space in 
variables (r,v±,(p,t), which we would still refer to as Q. 

Let o~(r, tp) be the distance of flight between two thermostats that originates at 
(r,<p). If, in addition, the flight originates with a given v±, then the speed during the 
flight is v j- ; . In order to show that this invariant measure is finite, we need to show 
that the time of flight °"( r '^) cos <» j s j n t e grable with respect to dpi, 



I v{r,<p) <**{*) f y) cos(^) ffa_v^_ , iV± ^) drdv ^ 

J v± J v± V vr cos 2 (v9j 

J V 7T COS{(p) 

<a ro ^ J f(r,v ± ) C os^(w))e-^> 2 dwd r dv ± 

//3max f f{r,V±) ft 

max / - 1 - , 

< o"max\/ / drdv± = a m! , x \ / — dp. 

V 7t J v ± V vr J v ± 

Here cr max is the maximum distance of flight and /3 max = max{/3i, • • • , f3 p }. 

The potential V introduced in [TT] was V(v±) = v^ 1 for v± < w™ in for some small 
v™ m . By Theorem 3.2 in [9] we know that J V(r, v±)dp < oo. We conclude that 
f ^-dp < f V(r, v±)dp < oo. Therefore, 
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/ diii < oo 

J v± 

and the invariant measure for the suspension flow is indeed finite. □ 
3.3. Mixing 

The proof of mixing and convergence of initial distributions to the invariant measure 
for $i relies on sampling the Markov process $ t at integer times, which generates a 
discrete-time skeleton chain $1. Denote the transition probability kernel for $x by V\. 
The following theorem by Meyn and Tweedie relates skeleton chains to the ergodicity 
of the Markov processes. 

Definition 3.1 The Markov process <& 4 is called ergodic if an invariant probability 
measure v exists and 

lim ||P*(z, •) - v\\ = 0, Vz e Q 

Theorem 3 fllfl Suppose the Markov process $t is irreducible and v is an invariant 
probability measure for & t - Then $ t is ergodic if and only if there exists a skeleton 
chain which is irreducible. 

The mixing of the invariant measure for the Markov process $ t and the convergence 
of initial distributions to the invariant measure follow from ergodicity by the Dominated 
Convergence Theorem. Indeed, 

lim - \i\ = lim sup | / (7>*(z, A) - fi(A))du\ 

< lim / \\P t (z,-)- f i(-)\\du = 0, 

and to show mixing one may replace Q with B in the second integral. 

Ergodicity follows once we show that the time-1 sampled chain $i is irreducible. 
Prop. 13.11 guarantees that for any two initial states (x,v) and (x',v') in Q, there exists 
a continuous time sample path from (x,v) to (x',v f ). 

We must now show that such a sample path from Prop. 13.11 can be modified to last 
integer time. Without loss of generality we can assume that each of the continuous-time 
paths from Prop. 13.11 passes through a given (a;o,fo). If there exists a family of closed 
sample paths from (x , v ) to (so, v ) such that the fractional parts of the times they 
use cover the interval [0, 1], then a particle can follow these sample paths to make the 
total time of traversing the path integer. By continuity, there exists a range of paths of 
this kind covering the time interval [ii,^] f° r some ti < t 2 . Then, if we repeat each of 
these paths k times, the time interval of the k th iteration covers [kti, ktz\. If k is large 
enough, kt 2 — kt± > 1, which means that all possible fractional parts of the times to 
complete such iterated paths cover [0, 1]. By adding an appropriate path of this kind 
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to the existing path from (x,v) to (x',v') through (x ,v ), we produce a sample path 
taking integer amount of time to complete. 

Similar to proof of Prop. 13.21 by modifying Props. 5 and 6 in [H] for any v±, v' ± while 
dropping lower bounds on the densities, we conclude that for all (x, v) and (x 1 , v'), there 
exists n and an open ball U( x >y) E5 (x f ,v f ) such that Vi((x, v), •) has strictly positive 
density on U^^y It follows that the time-1 chain is irreducible. □ 

3.4- Mixing is not exponential 

The main reason for sub-exponential rates of mixing rates is that, occasionally, the 
particle acquires low speed due to randomness of thermostats, which leads to long 
traveling times between collisions. We show that the (invariant) measure of the set of 
particles that would not collide with a thermostat before time r is of the order of 
This measure can be estimated using the information about the system used in the 
proof of the main theorem in [TT] . Then the dynamics is similar to the dynamics of an 
expanding map with a neutral fixed point (aka Pomeau Manneville) and can be treated 
similarly to [17J. 

Proposition 3.3 There exist (many) initial probability distributions A that converge to 
the unique invariant measure v with sub-exponential rates, i.e., 3q > such that for t 
large enough 

\\v:x-v\\>^ 2 . 

In particular, the unique invariant measure v for the Markov Process $ 4 is not 
exponentially mixing, i.e., there exist (many) Borel sets A, B such that 

| J[V T (z,A)-u(A)]dp\>u(A).^. 

B 

Proof of Prop. [3731 

Let B t be the set of all z e Q such that the particle will not collide with dT 
before time t. We would like to estimate v{B t ). In the equilibrium situation (when 
P = 1 = -.. = P P ), 

dv = c\ - % e-^ /cos2{lp) dv ± drd<pdt 

V 7T COS^ ((f) 

= cs 2 e~^ s2 cos(cp)dsdrd(pdt, 

where s = c ^ s ^ is the speed of the particle. The particle with speed s will not 
reach the thermostat in time r if t £ [0, CT ^ r '^ — r), so 



v{B T ) 




cs 2 e 



-/3s 2 



cos((p)dtdsdipdr 
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7r a(r,<p) 
2 r 




S 



ar -| o 



2 




ar-| 



2 




a(r,^) 0FrErf[^)] 
c — -= —3 J cos(<p)d<pdr 



[^^ + 0(\))cos(cp)d<pdr^\ 



9r 



where / ~ g means /(r) = 0(g(r)), i.e., there exists £, £' > such that 
?9(t) < f(r) < £g(r), for all r. 

We would like to show that the same result holds in the non-equilibrium setting. 

Lemma 3.1 There exist 5 > and tq large enough so that for r > tq, v{B t ) = 
v(B T+s ) + e(^). Moreover, v(B T ) « 

Assuming Lemma 13.11 we note that the dynamics for our system is statistically 
very similar to the dynamics of an expanding map with a neutral fixed point (aka the 
Pomeau-Manneville map). Indeed, the mass originally in B T evolves deterministically 
for time r; extra mass is 'deposited' from the parts of the phase space where particles 
experience nearly tangential collisions (which reduces the speed). We complete the proof 
of Prop. 13.31 using an argument very similar to [T7] . 

Let A <C v with dX = (pdv be such that ip > 1 + c on B To for some small c > 0. For 
t > r , let /„(t) = [v(B T ) - u(B T+nS )]/n. By Lemma [O, f n {r) rs ^ uniformly for all 
neN. i.e., ^ < f n (r) < £g . Then for k5 > r 



\V:\ - u\\ > (V:\)(B kS ) - u(B kS ) > X(B kS+nS ) - u(B. 



kS) 

> (1 + c)v{B kS+nS ) - u(B kS ) = [(1 + c) < Bk ^ _ l]u(B k s) 

is(B kS+ri5 ) + nf n (kd) 

= [c ~ (1 + c) , R n/ " ( f } f , M > (^) 
is(B k5+nS ) + nf n {kd) 

Then for A/" large enough and = iVn 

\\V* X - i/ > c - 1 + c — — — - ]v(B NnS ) > - =: — . 

v(B( N+X)nS ) +nf n (Nnd)) 2N 2 n 2 5 2 n 2 

□ 
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Proof of Lemma 13.11 In equilibrium, ((3 = Pi = ■ ■ ■ = /3 P ), dfx = 2v±e /3v ^dv±dr. 
Observe that the density of the steady state \x grows linearly near v± = 0. 

In non-equilibrium, we can obtain an upper bound on the density growth using the 
potential V in [llj. A simple computation ensures that the potential V in [TT] can be 
modified to be V(r,v±) = 1/fJ for v± < t>™ m and any < 7 < 2. By Theorem 3.2 
in [9] , fV(r,v±)dfjL < 00. Therefore the density of /x must grow at least linearly near 
v± = 0, i.e., f Qr f(r, v±)dr < ('v± for some 

Let us obtain an upper bound on u(B T ). The maximal particle flight time is a max , 
so the maximal possible speed in B T is £m ^ > v j r . i.e. v± < £saz - 



CT max 

2 CTmax 



v{B T ) < J J J J cf(r,v ± )p v± (ip)drdtdpdv ± 



-| ar 

<y max tt 

T 2 



<-u,ax / ! cCvJ^-^e-^^d V dv ± 
V 7T COS^(</9) 



< If / / v 2 ± e-^ v > 2 dwdv ± = K I v ± dv ± < — 



2 r 2 r 2 



-00 
We will obtain the lower bound using the following lemma 



Lemma 3.2 Let /ii 6e the lift of fi to the (r,v±,ip)- space and let A(v±) := {(r, v±,(p) : 
^mm < v± < - u ™ ax ) -y^ < T/ien there exists v° > suc/i t/iat 

for vj_ < v° ± and some K > 0. 

We will prove Lemma 13.21 at the end of this subsection. 

Note that if ~v\ <C t>™ in and (r,v±,ip) G then is uniformly bounded 

above by some $0 > and the time of flight between thermostats is bounded below by 
a min cos(^ )/^r x - 

Let A(v±) = {(r,v±,ip,t) G Q : (r, v±,(p) G A(v±)} and fix 5 < a min cos(i9 )/^j nax - 
Assume at time a particle originates at (r,v±,(p,t) G and is going to collide 

with a thermostat within time 5. Assume that at collision, the new angle ip new drawn 
from p v > ± ((Pnew)d<Pnew satisfies |y nem | < f . Then the particle's speed after the collision 

ls s ' = cJat-) < ^- L ' ^ ^-L < 2(t+5) anc ^ Wnew\ < f, then, after the first collision, the 
particle will not collide with the next thermostat in time r + 5 and the coordinates of 
the particle belong to B T at time S. 

Then, by Lemma 13.21 the //-measure of states that get to B T in time 5 is greater or 
equal to 

7T — £ 

W(\<Pnew\ < o ) X U ( A (V±)) > C5KV\ =1 -\ 
3 " T 6 
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Thus v(B T ) > v{B T+5 ) + £ and u(B T ) > £ 77 ^ F « £ □ 
Proof of Lemma 13.21 

Fix and t^ ax = Let C := {(r, vjj : ^ in < v± < w^ ax }, nad let m c be the 
uniform probability measure on C. The minorization condition [1] states that: 3N and 
rj such that 

inf V N ((r,v±), ■) > rjmc- 

{r,v ± )eC 

Therefore, by invariance, \x has density at least fJ*(C)r) on C. Note that /x(C) > since 
the chain is irreducible. 

Given (r,v±), we would like to estimate the probability that v' x = ° c ° s s ^ v± < v±. 

Using the same parametrization as in subsubsection 4.2.1 in [TT] we find that ^° s s2 ^ is a 
piecewise quadratic function when viewed in w = tan(<^)-coordinates. We are interested 
in integrating the angle density p v± (ip) on the intervals (w ,Wi) := (tan(</? )> tan(<^i)) 
on which < < Sl. 



P« < «i) = / J^^—e-^^d^ = / J^v ± e-^> 2 dw 

J \ 71 COS Z (ip) J V 7T 

cosfi/j') ^ nj_ cos(i/)') ^ 

cos((/j) — 1>J_ cos(i/?) — Vj_ 

= S(Erf[v^t; ±Wl ] - Erf[VA^o])/2. 

The difference in the error functions is only substantial when y/j3~iV±wi and yf]3iV±wo 
are sufficiently close to 0. In particular, this happens near geometric locations r e dT 
from which a radially emitted particle collides with the next thermostat tangentially. 
See Fig. [2j We will show that, for small enough v±, the measure of {r : such that 
wo < < wi (or wi < < wo)} is > t(y±) 2 for some t > 0. 

Let ro and r be such that if a particle originates radially from ro (with (fo = 0), 
then it arrives to r tangentially (i.e. (p' = ±| ). Note that <^o = corresponds to 
ifo = 0. Now we can start shifting r keeping w < < wi (or ui\ < < w ) until we 



Ar 

= u, i.e. cvayyji) = • J-ue leiciiiuii ueiween '</v = 

follows: 

i?i — (To tan('0) i?i 
sin(( y 9 / 1 ) cos(?/>) 

Solving for sin (ip) we get 



get W\ = 0, i.e. cos(v9 / 1 ) = ^kfe. The relation between tj) = ^ (see Fig. and ^ is as 



-Ri cos(^) — <t sin(^) = .Ri sin(<^ 



-o-qRi sin(v?' 1 ) + R l cr Jl + -4 cos 2 (v?i) 
Bin(^) = -r-f 2 ^ 



-a 0J R u /l - (^W) 2 + ^oa/I + §1 
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Figure 2. Geometric locations that produce small v\ 



> „2 i d2 a (tss) = :i K) 



From computations in subsubsection 4.2.1 in [TT] one obtains that 2 ^~ x — l w i ~ 
u^ol < Therefore, for Uj_ is small enough, if r is chosen such that u> < < u>i (or 

wi < < X w ) and < u ± < ^ ax , 



P« < t?±) > Erf[VAu±wi] - Erf[VA«±«'o])/2 > Erf[VA«i>i - ™o|]/2 > , 

V 7T 4(To 

Therefore, ^(A(v ± )) > ^C) V l(v ± ) 2 V ^^ =: ^ x . □ 
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